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1. INTRODUCTION 

Brushless direct current (BLDC) motors are one of the popular applications in industrial sector. In 
particular, this motor is used to control the servo motor systems and its electrical equipment because of its 
robustness, high efficiency, produces high torque for volume capacity ratio (VCR), and simplicity [1], [2]. 
The BLDC motor is studied for diverse energy levels at diverse applications. In order to get the parameters in 
BLDC motor, deploy Computer Aided tools but this tool takes long time to get this parameter in BLDC 
motor. Hence, in this article, we have introduced an Evaluation and tested system model to get the original 
parameters while designing BLDC motor. To optimize original parameters in BLDC motor, use MAXWELL 
3-D FEM (finite element method) [3]. 

The effect of electrical motor is not working with model test continuously since it needs high 
consistency, low cost, small in size, high pull-out torque to load and to activate in manufacture in industrial 
sector [4]. In this article, the angle which is present in the lower torque is chosen by rated speed and its speed 
while designing the motor in optimized manner. Hence, this optimized design is used to enhance the potential 
of overload while handling this design efficiently and get large torque during short time. Hence, this produces 
the BLDC motor design in regular manner and also reduces design time [5]. In the account of design 
consideration, the initial design parameters are magnet assembly gap, rotor saturation, Assumed flux leakage 
coefficient, rotor saturation, and permanent magnet [6]. 

Interior permanent magnet brushless DC (IPMBLDC) electric motor is the main type of these 
motors, constructed with the permanent magnets inserted into the rotor core [7]. The pathway of leakage in 
an interior permanent magnet motor comprises a saturable magnetic reaction which creates the coefficient 
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variable for flux leakage [8]. This article determines and to test the design in the interior portion of 
permanent magnet bruless DC motor and it is used to concentrate the windings and also, apply in electric 
parameters in this motor. The equivalent magnetic circuit of the IPMBLDC motor is to analyze the 
coefficient of flux and working selected point of a permanent magnet based on the design output parameters. 
The MAXWELL 3D FEM is deployed to validate the equivalent circuit of magnetic field and optimize the 
IPMBLDC motor design. To validate the design, the experimental responses are shown. 


2. MOTOR DESIGN PROCEDURE 

The interior permanent magnet BLDC motor design procedure is very difficult when compare with 
other motor [9]. In IPMBLDC motor, the leakage which is produce in magnetic flux and its path generally 
comprises saturable magnetic circuits. The BLDC motor of 24 slot IPMBLDC design procedure as sown in 
Figure 1. The flux leakage coefficient ọL and the operational point in the permanent magnet wm are 
deployed for the first design process in IPMBLDC motor. The metric value for ọL and Ym are assumed after 
selected parameters are evaluated roughly based on our equivalent circuits. We offer the primary design size 
of the IPMBLDC motor, system model, and a small torque angle value value when compare to legacy 
designs by the rated load. Therefore, we get greater overloading handling capacity and enhanced total 
efficiency. Then, to verify the design process and parameters, analyse the finite element and simulate FEM. 
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Figure 1. 24 slot IPMBLDC design procedure 
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In Figure 2, shows the configuration of 24 slot IPMBLDC motor. The loop of core flux flows via the 
magnet, yoke in stator, air-gap, winding in stator, tooth in stator, and yoke in rotor. The flux linkages taking 
in to consideration for permanent magnet, the IPMBLDC motor equivalent magnetic circuits and these two 
combined with half-pole pair to assume the equilibrium for this motor. 


Rs 
R = ar Ret + Rag + Rra + Rro (1) 


Where Rl is the total reluctance of air gap; Rs is the reluctance of stator yoke; R,; is the reluctance of stator 
tooth; Rag is the reluctance of Air gap; R,q is the rotor reluctance above the magnet; R,, is the rotor 
reluctance below the magnet. 
The magnet ending flux leakage reluctances can be expressed as (2). 
dg 


Ro = (2) 


HoAg 


Where do is the magnet distance and the yoke channel Ac is the area of the cross section in the air gap. The 
magnetic reluctance can be expressed as (3). 
— _Hm 

HoHrAm 


Rmo (3) 
Where hm is the permanent magnet length. Am is the area of the cross section in the magnet. The summation 
of magnet ending flux leakage reluctances and reluctance of magnet can be expressed as (4). 


Ro = Ro + Rmo (4) 


To calculate the consistency of the certain operation point and the coefficient of the flux leakage, it is 
essential to analyze the equivalent of the magnet circuit [10] as shown in Figure 3, we consider the density of 
web flux in the magnet is 1.8 tesla. Hence, the leakage of the flux in magnet ọL is obtained. The density of 
flux in the magnet intensity is limited to 2 Tesla. 


ory 


Permanent Magent 


Figure 2. 24 slot IPMBLDC motor configuration Figure 3. Magnetic equivalent circuit 


3. THE DESIGN AND CALCULATION OF BLDC MOTOR 
3.1. Dimension of the permanent magnet 

When the air-gap dimension is obtained, the (5) and (6) can be used based on [10] to conclude the 
dimension of the permanent magnet. Assumed flux leakage coefficient and selected permanent magnet 
operating point are used in the equations. 


— KsKx%m 


~ (1m) PL 


m (5) 
Where hn is the permanent magnet length; g is the air gap length; K, is the motor saturation factor; Kx is the 
rotor structure factor. The IPMBLDC motor air gap flux produces square waveform, hence the metric rate of 
the pole is coefficient must be greater for sufficient case see in [11]. The permanent magnetic width can be 
expressed as (6). 
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Where Bg, is the peak value of air gap in fundamental wave; Kọ is the air gap flux waveform factor. In 
IPMBLDC motor, the coefficient for pole is formulated in (6). The metric value of the saturation factor in 
motor varies between 1 and 1.3. The value of the structure factor for rotor between 0.7 and 1. 

The permanent magnet material B-H curve is shown in Figure 4(a) (see Appendix). The permanent 
magnet operational point is related with motor energy. Bm is the flux density and Hm is the magnetizing force 
in the magnet. There are two lines in the operational point of the permanent magnet known as load and 
device line. The line in the load is mentioning the line via the starting point to the operational point and the 
device line is distinct as an ideal demagnetization curve [12]. The ideal IPMBLDC motor air gap flux 
waveform is square. Some derivation, the permanent magnet flux density is formulated as: 


Hm= - (1- wm) He (7) 
Bm = -(uoHm) (8) 
Pc = (Ag/g)(Am/dm) (9) 


Where Pc is permanence coefficient; Ag is air-gap area; g is air-gap length; Am is cross-sectional area; dm is 
the distance between the magnet; wm is called the operating point of a permanent magnet. 

The simulated B-H curve of the permanent magnet material is shown in Figure 4(b) (see Appendix). 
In this experiment we have choosen carbon steel. It is a normal carbon steel with 0.10 percentage of carbon 
content. We achieve the selection of the power in maximum level is ~wm=0.5. In order to overcome the 
demagnetization effect, we always select an operating point larger than 0.5. In the proposal, the operating 
point between 0.5 and 1 is chosen [13]. To find out the air-gap of the magnetic field from the rotor of PM. 
We have used as an expansion of fourier series; the air-gap of flux density in magnetic field is given by: 


> | 22M n+2PM 
Bi= = Sve Bm cos(he) + J,-Pbu —B,n cos(he) (10) 
5 es a1) 
Kpn = sin(h * (12) 
By pe= = sin CE) By (13) 


Where Kpn is the pitch factor of hth harmonic; B, is the air-gap flux density; B} px is the peak air-gap flux 
density. 


4. PERFORMANCE INVESTIGATION 

To validate the design parameters and magnetic circuit model, the 3D Finite element investigation 
has been used. We consider the density of the flux is constrained from 1.6 and 1.8 Tesla. For this IPMBLDC 
motor, Magnet 3D is deployed to estimate and improve the original factors are obtained from the equation for 
motor design. The design purpose is to decrease the ripples which are present in the torque and enhance the 
efficiency. Certain factors namely air-gap length, magnet length, and so on will be a little modified to get 
improved performance [14]. The dimensions of the IPMBLDC motor model are shown in Table 1. 


Table 1. Dimension of the motor model 


S.L. Model in motor Range 
1 The slots/pole Numbers 12/4 
2 Outer diameter/mm in Stator 48 

3 Inner diameter/mm in Stator 24 

4 length/mm for Air-gap 0.2 

5 Length/mm in Magnet 9 

6 Voltage in volts 200 
7 Current in Amps 1 

8 Power in watts 200 
9 Speed in rpm 2500 
10 Torque in N-m 1 
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In Pulsating, torque is the relation impacts between the permanent magnet field in the rotor and the 
teeth of the stator pole [15]. It builds a torque pulsating that does not contribute to the effective of the total 
torque [16]. The flux distribution waveform angle and torque characteristics for the IPMBLDC are shown in 
Figure 5 and Figure 6. The ampere-flux distribution of the stator residue is constant and predetermined in 
space for a fixed commutation interval as the magnet rotates the previous period, produce a linear phase 
difference variation of flux-linkages, and from flat-top back-waveform for EMF [17], [18]. To get high 
torque density, modification done in stator concentrated windings and rotor dimensions are used [19], so the 
value of torque ripples and motor efficiency is suitable for 12-slots, 4-pole motor. The speed characteristics 
of the IPMBLDC motor model are shown in Figure 7. The speed versus torque plot of the IPMBLDC motor 
model is shown in Figure 8. The torque increases, the motor speed is decreased. The plot shows that the 
design model has around close to the theoretical existing motor characteristics graph. Detailed simulated 
relationship between torque, flux density, and speed is shown in Table 2. 
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Figure 7. Speed characteristics 
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Table 2. Detailed simulated relationship between torque, flux density, and speed are shown 


1090 


1500 
Speed in RPM 


2000 


Figure 8. Speed vs torque characteristics 


2590 


Time in sec. Flux Linkage Speed(deg/sec) _ Load Force/Torque(Nm) 
0 -0.008011516828 0 -0.6 

5 -0.006536902708 -710.6986868 —-0.6303898785 
10 0.0005355561372 -1242.931520  -0.6531484561 
15 0.003531455894 -1607.503179  -0.668737746 
20 0.005934708220 -1822.05293 -0.6779120158 
25 0.000345140740 -1966.694286 ~— -0.6840969621 
30 -0.006862879536 -2078.011213 -0.6888569370 
35 -0.002448606938 -2151.676149 -0.6920068914 
40 0.005540622231 -2215.866840 -0.6947517217 
45 0.002907926839 -2261.638092 -0.6967089264 
50 -0.005417490002 -2308.086414  -0.6986950830 
55 -0.002690250688 -2343.616719  -0.7002143789 
60 0.005971689114 -2376.527811 -0.7016216759 
65 0.00190826119 -2395.715475  -0.7024421513 
70 -0.005346391726 -2414.224014  -0.7032335869 
75 -0.000732794433 -2434.628720 — -0.704.1061037 
80 0.006230476736 -2443.034874 —-0.7044655556 
85 -0.000324799638 -2456.101107 -0.7050242751 
90 -0.006544184562 -2464.358617  -0.7053773709 
95 0.002230922776 -2475.212666 —-0.7058414962 
100 0.00592900262 -2488.032278  -0.7063896701 
105 -0.003704002075 -2485.111041 -0.7062647297 
110 -0.003651447474 -2483.060521 -0.7061770750 
115 0.00543777656 -248 1.82806 -0.7061243745 
120 0.001876431442 -2479.914360  -0.7060425433 
125 -0.005751353033 -2480.390204  -0.7060628907 
130 -0.000106910481 -2483.845997  -0.7062106623 
135 0.006525068983 -2484.412592  -0.7062348903 
140 -0.001449741492 -2488.078118 -0.7063916302 
145 -0.006225249296 -2493.286151 -0.7066143287 
150 0.003320757915 -2495.135609  -0.7066934126 
155 0.003953507962 -2499.779089  -0.7068919704 
160 -0.005128059792 -2504.419074  -0.7070903787 
165 -0.002124128813 -2505.476122 -0.7071355787 
170 0.005605705254 -2507.624183 -0.7072274311 


3000 


161 


The analysis of the fast fourier transform (FFT) in flux distribution of air gap is shown in Figure 9. 


Only odd harmonics ripples are appeared but even harmonics ripples are disappeared, which results that the 
ripples of flux density distribution in the air-gap are in excellent condition [20]. Measuring the torque of the 
motor is very difficult. It requires costly equipment. So we suggest calculating the efficiency of the motor. It 
is calculated as the ratio of mechanical output power to electrical input power. 


E = Pout 


Pin 
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Pout = Pin XE (15) 


Output mechanical power of the motor is calculated by: 


Pour =T Xw (16) 
@=rpm x = (17) 
Pous =T X rpm x = (18) 
2T 
Pn XE=T Xrpm X (19) 
IxVxE=T xrpm x = (20) 
T= IXVXEx60 (21) 
rpm x 2T 


Efficiency is chosen between zero to maximum; in our simulation below 2500 rpm is optimal speed for 
calculations, it is used 100% efficiency (E = 1). 


I =1A4,V = 200V,E = 1 and rpm = 2500 
T = (1 x 200 x 1 x 60)/(2500 x 2 x 3.14) = 0.76 N -m 


Where Pin is the input power in watts (W); I is the current measured in ampere (A); V is the voltage 
measured in volts (V); T is the Torque in Newton-meter (N-m); is the angular speed in radian per second. 
The torque efficiency is based on equation (21). For the simulation efficiency is 93%, some mechanical 
losses are neglected and the parameter is chosen the very nearest and rounded value, this is the reason 
simulation efficiency is 7% lower than the theoretical efficiency. 
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Figure 9. Fast fourier transform analysis 


5. CONCLUSION 

In this paper, we compared to design the IPMBLDC motor using 3D-Finite element analysis and the 
structure of the theoretical design for electrical effect determination of IPMBLDC. The representation of the 
magnetic circuit has been developed to evaluate the permanent magnet dimension, the operational point of a 
permanent magnet, coefficient of flux linkage, and torque efficiency. The theoretical design has been 
compared and confirmed using MAXWELL 3D analysis based on FEM. Finally, torque efficiency is 
compared with the theoretical design structure and analytical method of FEM. The theoretical calculated 
torque density is 0.76(N-m) and simulated torque density is 0.70(N-m). The simulated torque efficiency of 
93% is very close to the theoretical torque efficiency. Therefore, the MAXWELL 3D analysis tool is most 
suitable for optimized and analyzing structure and evaluation of electrical impact 
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Figure 4. The permanent magnet material: (a) B-H curve of the permanent magnet material and 
(b) simulated B-H curve of the permanent magnet material 
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